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Abstract: 

In this paper we consider high-intensity source of electron antineutrinos from the production and subsequent decay 
of 8Li. It opens a wide range of possible searches for beyond standard model physics via studies of the inverse 
beta decay interaction. In IsoDAR experiments Lithium 8 is a short lived beta emitter producing a high intensity 
anti-neutrinos, which is very suitable for making several important neutrino experiments. In this paper we used the 
GEANT4 program, to simulate neutrino production using proton and deuteron beams. We find that the neutrino 
production rate is about 3 times from deuteron beam than from proton beam in low energy region. 
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1 Introduction 


2 IsoDAR Design 


As the next step, neutrino physics requires high in¬ 
tensity experiments for precision measurements. We pro¬ 
pose the neutrino experiment on proton or deuteron 
beam hits target. The Isotope Decay At Rest (Iso- 
DAR) n 0 technique provides a high-intensity, low en¬ 
ergy source of antineutrinos with sensitivity to antineu¬ 
trino oscillations. The experiment can perform com¬ 
pelling tests of models for new physics that explain high 
Am^ oscillations through the introduction of one or more 
sterile neutrinos. 

IsoDAR can emit anti-electron neutrino at a rather 
high-flux by beta decay of short lived isotope. Lithium 
8 is a nice candidate of this isotope and it can be pro¬ 
duced by injection of proton or deuteron on a beryllium 9 
target. To meet the required neutrino flux of oscillation 
experiment, the intensity of proton beam should be at 
least 10 mA with energy of 60 ^ 1000 MeV. This require¬ 
ment is the main limitation to IsoDAR source for the last 
decades. But now, by taking advantage of high intensity 
proton beam technique, this becomes very promising and 
worth a detailed simulation. IsoDAR sources provide an 
opportunity for the precision measurements required for 
the future of neutrino physics [3]. 



Fig. 1. Initial design for IsoDAR target. 


Fig. [T] shows initial design for IsoDAR target 
(Ref. [4]). It’s like a neutrino converter: absorb proton 
(or deuteron) and then emit neutrino. 

When proton inject on the Beryllium target, several 
interactions happen subsequently: 


p + ^Be — ^Li^2p 
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3 Simulation Results 

p^^Be — >n-\-p-\-2 ^He (2) In the present work, the GEANT4 simulation code 

[3] has been used to simulate the low-energy protons in¬ 
duced isotopes production in Be targets. 

After a simulation of a number of proton incidence, 
Li —Li (3) we can get the neutron spectra crossing different barrier 

components as Fig. shows: 


^ Li—y v^-\-Be (4) 

With proton beam incidences, the Beryllium 9 tar¬ 
get can produce Lithium 8 and neutron, to which we 
pay more attention, as Eq. [^and Eq. describes. Also, 
we can see from Eq. that the production of neutrino 
is directly related to that of Lithium 8. So in our pri¬ 
mary design, there is a FLiBe sleeve outside and inside of 
the Be target respectively to absorb neutrons and pro¬ 
duce Lithium 8, which process is described by Eq. 
FLiBe is molten salt made from a mixture of lithium 
fluoride (LiF) and beryllium fluoride (BeF 2 ) with abun¬ 
dance of Li 7 up to 99.99%. It has a high melting point 
and a better mechanical property compare to elemental 
Lithium. The heave water between Be target and inner 
FLiBe sleeve is moderator for neutron to improve its ab¬ 
sorptivity by Lir, besides, the heavy water sleeve can be 
design as cooling system to carry out heat from beam 
injection. The graphite layer outside the outer FLiBe 
sleeve is set as a reflector to improve the utilization for 
neutron. To prevent neutron radiation, a thick concrete 
sleeve is set outside covering the whole equipment as a 
shielding layer. 

While using Geant4 [6] to simulate interactions of 
neutrino production, we need to specify relevant model 
for every process and every particle, i.e., the PhysiesList 
package in Geant4. In our primary simulation, we choose 
QGSP_BIC_HP, a reference physics list from Geant4, 
to describe relevant reaction channels and cross section. 
This package implies quark gluon string model to high 
energy interaction (5-25 GeV), the excited nucleus cre¬ 
ated is then passed to the pre-computed model modelling 
the nuclear de-excitation. As for interactions below 10 
GeV, it uses binary cascade model (BIG) to describe the 
propagation of incident hadron and secondary through 
nucleus, which is regard as a series of two-particle colli¬ 
sion. For energy from 20 MeV down to thermal particle, 
this model takes data driven high precision neutron pack¬ 
age (NeutronHP) to transport neutrons below 20 MeV 
down to thermal energy. Gompare to other physics list 
package, QGSP_BIG_HP is able to rebuild the produc¬ 
tion of secondary particles from p—A or n—A interactions 
well and is more precise in low energy region. It can gives 
us a more accurate isotope production rate and neutrino 
distribution in IsoDAR. 


Neutron spectra 



Fig. 2. Spectra for neutrons crossing the barrier 
components when 10^ protons incident with en¬ 
ergy of 60 MeV. 


From the above figure we can see the following points: 

1. A mass of neutrons are produced after protons inci¬ 
dence. And the neutrons coming from Be target to 
heavy water sleeve has the highest average energy. 

2. We can see the moderator effect of heavy water. 

3. Many neutrons escape out from the target sleeve, 
so it’s necessary to add a thick FLiBe sleeve sur¬ 
rounding target to have a good usage of neutrons. 

Besides, we can get the isotope production in differ¬ 
ent parts from proton incidence and deuteron incidence. 


inner FLiBe sleeve 
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Fig. 3. Isotope production in different parts. Top: 
60 MeV proton incidence; Bottom: 80 MeV 
deuteron incidence. 


Fig. a shows the production rate of Neutron and Li8 
respectively with different incident energy for proton and 
deuteron beams, we can see approximately a linear rela¬ 
tion between production rates and incident energies, We 
fit the two plots and get the following formulas. For the 
neutron production rate, 

A^P^ = 2.7xlO-^xF;^•^ (5) 


A^p, = 1.1x10-"xF;'-\ (6) 

and for Li8, 


We can see from Fig.j^that Li8 is mainly produced in the 
outer FLiBe sleeve, which is reasonable as more neutrons 
enter into this sleeve. Also, we can see that Deuteron’s 
Li8 production rate is higher than Proton’s. This result 
is meaningful and worthy of further research as it will 
tell us which incident beam is better in IsoDAR experi¬ 
ment. So we have a comprehensive comparison between 
proton and deuteron incidence with different energies. 


Neutron production varying with incident energy 




Fig. 4. Neutron (Up) and Li8 (bottom) produc¬ 
tion rate at different incident energies, for proton 
and deuteron beams respectively. 


7Vpp,8 = 1.4xlO-"F;'•^ (7) 

A^i.Lz8=9.1xlO-"xF;°-^N (8) 

Here, Np^ and corresponds to the neutron produc¬ 
tion rate for proton incidence and deuteron incidence 
respectively, while Nppis and Nppis means the Li8 pro¬ 
duction rate. R^s in the above plots are close to 1, which 
means this linear fitting is pretty accurate. Compare pro¬ 
ton’s production rate at 100 MeV and deuteron’s produc¬ 
tion rate at 200 MeV, we can see that the later is not 
simply twice as many as the former. This indicates that 
our simulation includes the nucleus effect, which is more 
realistic. 

To have an intuitive comparison between proton and 
deuteron beams, we have this picture below. 


03 
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Fig. 5. Li8 production ratio (Blue) and neutron 
production ratio (Red) of Deuteron injection over 
Proton injection. This two sets of data have been 
fitted with a power function: ratio = po x 


In Fig. we compare the neutron and Li8 production 
rate for proton and deuteron incidence with different en¬ 
ergy, and then fitted them with power function respec¬ 
tively. We can see that in the low energy region around 
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100 MeV, the Li8 production rate for deuteron is about 
three times of that for proton. But as beam energy in¬ 
creases, the ratio approaches to about 1.5. So we can 
expecet that several hundred MeV deuteron beam has 
much advantages than proton beam with half energy to 
generate neutrino. 
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Fig. 6. The 8Li isotope DAR anti-electron- 
neutrino flux spectrum distribution. 


vV+v 
+ + + + 


Entries 100000 

Mean 6.72 

RMS _ 2.533 


We simulate one hundred thousand Li8 decay at rest. 
The neutrino spectrum produced from IsoDAR Li8 is 
shown in Fig. we can see that the spectrum is in 
agreement with experimental result [7]. 


4 Discussions and Conclusions 

IsoDAR represents both a novel concept of applica¬ 
tion to neutrino physics measurements. Neutrino exper¬ 
iment on the China Accelerator Driven Sub-critical Sys¬ 
tem (ADS) is being proposed [8]. The first phase of ADS 
is named China Initial ADS (CIADS). The CIADS facil¬ 
ity will offer the highest intensity protons or deuteron 
beams over the world with energy between dozens of 
MeV and 1 GeV. 

IsoDAR is a novel design, extensive and precise sim¬ 
ulations targeting the most challenging aspects of ex¬ 
periments are highly required. In this preliminary sim¬ 
ulation, we can see that by taking advantage of ADS 
technique, IsoDAR is a feasible neutrino source for ex¬ 
periment searching for sterile neutrino. And we recom¬ 
mend deuteron beam to hit the target as it has advan¬ 
tages over proton below 1 GeV. As we see, the deuteron 
beam has much advantage over proton beam around sev¬ 
eral hundred MeV. It will offer an opportunities for high 
precision neutrino measurements, such as sterile Neu¬ 
trino Searches, precision electroweak tests of the stan¬ 
dard model, and coherent neutrino scattering, etc. 
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